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1. Introduction 

Many authors emphasized difficulties encountered in the separation of the total angular 
momentum of light into its orbital and spin parts [Tj-[18j. A popular formula 
expressing this separation, presented in several textbooks [H]-[2T] and used in pQ-[5], 
has the form 

Jd 3 rrx (e E x B) = J d 3 re E 2 (r x V)^ + J d 3 re E x A. (1) 

This prescription is marred by a defect: the splitting is gauge dependent because 
it involves the vector potential A. This problem has been resolved by an ad hoc 
postulate that the potential must be evaluated in the transverse gauge but this 
prescription lacks a deeper foundation. In other papers the separation has been 
given only for monochromatic fields or in the paraxial approximation. An additional 
problem that has not been resolved to the satisfaction of many authors was caused 
by their wish to disentangle completely the orbital and spin degrees of freedom. This 
is possible for massive particles but not for massless particles. The direction of the 
spin for all massless particles is firmly locked onto the direction of momentum: it 
can only be parallel or antiparallel to momentum. In other words, the helicity of 
massless particles — the projection of its total angular momentum on the direction of 
momentum — can only take the values ±s. This fact makes it impossible to rotate 
independently the orbital and spin degrees of freedom of photons. 
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The correct, gauge invariant separation of the total angular momentum into its 
orbital and spin parts has been proposed long time ago by Darwin [32| . The Darwin's 
classic paper was cited without a comment only in [31 [S]. In several papers from 
our list [5], [H]-[IZ] the authors rederived his result, usually in a special case of 
monochromatic waves. The Darwin formula is based on the Fourier transforms of 
the electromagnetic field. Therefore, it does not suffer from gauge dependence. With 
slight changes of notation it reads: 

Jd 3 rr x (eo-E x B) = -2ie [E*(k)(k x V k )E t {k) + E* {k) x E{k)]{2) 

where E(k) is the plane-wave component of the electric field, 

/d 3 k 
^ m [E(k)e- i ^ +ik -r+c.c.}, (3) 

In the present paper we follow in the footsteps of Darwin who wrote "The main 
principle of the present work is the idea that, since matter and light both possess 
the dual characters of particle and wave, a similar mathematical treatment should be 
applied to both, and that this has not been yet done as fully as should be possible". We 
show that, indeed, the wave-particle duality enables one to determine the correct 
separation of total angular momentum. Namely, we shall show that the Darwin 
separation of the total angular momentum for an arbitrary electromagnetic field into 
two parts follows from the photon picture of the electromagnetic field. It is in essence 
the separation into the part perpendicular to the photon momentum and the part 
parallel to the photon momentum. The first part must be identified with the orbital 
angular momentum whereas the second part must be identified with spin — represented 
by helicity. In this way, by seamlessly joining the particle and the field aspect of 
electromagnetism, we complete the program started by Darwin. Our analysis of 
the angular momentum of light starts from the quantum mechanical description of 
photons. 

2. Quantum mechanics of photons 

There is no consensus as to what represents the photon wave function in the coordinate 
representation (cf., [33]) but there is no disagreement as to the meaning of the photon 
wave function in momentum space. This wave function has been introduced in the 
early years of quantum electrodynamics by Fock |24j and was used as a standard 
concept in textbooks and monographs [25 -[29 . Once we accept the existence of 
the photon wave function in momentum space we should define the action of various 
operators representing physical quantities. 

In a relativistic theory — and there is no nonrelativistic theory of photons — we 
should first of all define the operators representing ten generators of the Poincare 
group: the generators of translation in space (momentum P), translation in time 
(energy H), rotation (angular momentum J), and Lorentz boosts (moment of energy 
K). These operators must obey the following commutation relations appropriate for 
the Poincare group: 

[H,Pi]=0, [H,Ji} = 0, [H,K i ] = -ihcP i , (4a) 
[Pi, Pj] = 0, [J h Jj] = ihe ijk J k , [Ki,Kj] = -itk?e ijk J k , (46) 
[Ji,Pj] = ihtijkPk, [Ji,Kj}=iHeij k K k , [K t , Pj] = ihd tJ H . (4c) 
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There are no problems with the construction of the generators for massive particles. 
The following set of operators was given long time ago by Foldy [30] 

H = E p , (5a) 

P = p, (56) 

J = ihV p xp + S } (5c) 

^ = ihE ^-^TE- p <™ 

where V p denotes the gradient with respect to the components of momentum and the 
spin vector S is built from three (2s + l) x (2s+l) matrices that obey the commutation 
relations of angular momentum. The matrices S act on the (2s + l)-component wave 
functions describing the states of a particle with spin s. In this case, the splitting of 
the angular momentum into its orbital and spin parts is quite obvious. 

The representation of the generators of the Lorentz group for massless particles 
was given by Lomont and Moses [3T]. We will use here a modified version of these 
generators [571131] that exhibits its geometrical meaning. The momentum operator, by 
definition, acts on the wave functions in momentum representation as a multiplication 
by hk. There is no question that the operator representing the energy of the photon 
(the Hamiltonian) must be the modulus of the momentum vector multiplied by c. The 
complete list of generators contains also the operator of angular momentum and the 
boost operator, 

H = tkJk, (6a) 

P = hk, (66) 

J = \hD x k + hxrik, (6 c) 

K = ihuj k D, (6d) 

where n*. = k/\k\, the photon helicity operator x has the eigenvalues ±1, and D 
stands for the covariant derivative on the light cone (V*. = d/dk), 

D = V k -i X a(k). (7) 

These operators act on the two-component photon wave functions 

«*> - ( f%\ ) • < 8 > 

and satisfy the commutation relations (|4al) - (|4cp appropriate for the Poincare group. 
The two components of the photon wave function correspond to two eigenvalues of x, 

h(k) \ ( h(k) \ ( ) 

X { fn(k) ) { -fn(k) ) ■ (9) 

We used the indices L and R to denote the eigenfunctions of the helicity operator since 
they correspond to left-handed and right-handed circular polarization. The properties 
of the covariant derivative are obtained from the commutation relations for the angular 
momentum and they read: 

[D i ,D j ]=i X e ijl n l /\k\ 2 . (10) 

These conditions determine the vector ct(k) up to a gauge transformation 

a(k)^a(k) + V k v(k), (11) 
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which is connected to the change of the phase of the wave function, in analogy to the 
theory of charged particles coupled to electromagnetic field. The generators (pat - ([6d|) 
are hermitian with respect to the following Lorentz-invariant scalar product 

/d 3 k f d 3 k 

_ft(fe). fl (fe) = J —[fi(k)g L (k) + r R (k)g R {k)}. (12) 

3. Electromagnetic field 

In order to solve the problem of the total angular momentum separation into two parts 
for the classical electromagnetic field, we shall employ the correspondence between 
the fundamental physical quantities (energy, momentum, and angular momentum) in 
photon quantum mechanics and in Maxwell theory. In quantum mechanics of photons 
these quantities are represented by the operators (IBa|) - (|6dj) . In Maxwell theory these 
quantities are given as space integrals of corresponding densities built from quadratic 
expressions in field vectors. A very convenient tool in this construction is a complex 
vector F, 

F = ] l^(E + icB), (13) 

that was named the Riemann-Silberstein (RS) vector in [23]. The Maxwell equations 
expressed in terms of F are: 

d t F(r,t) = -icV x F(r,t), V-F(r,t) = 0. (14) 

The field energy H, the field momentum P, the field angular momentum J, and 
the field moment of energy K can all be constructed from the energy-momentum 
tensor of the electromagnetic field. These quantities expressed in terms of the RS 
vector are: 

H=^Jd 3 r[e a E 2 + B 2 /^] = Jd 3 rF*-F, (15a) 
P = Jd 3 r [e E x B] = ~ J d 3 rF* x F, (15b) 
J = J d 3 r r x [e E(r) x B(r)} = A J d 3 r r x (F* x F) , (15c) 
K= \Jd 3 rr [e E 2 + B 2 /^ Q ] = Jd 3 rr(F*-F). (15d) 



These quantities, like their counterparts in photon quantum mechanics (|6a[l - (|6rijl serve 
as the generators of Poincare transformations of the electromagnetic field. They 
have analogous algebraic properties of the Poincare group (|4ap - (|4cl) . with quantum 
commutators replaced by Poisson brackets, [a, b]/\h — > {a, b}, (cf., for example |27j). 

All solutions of Maxwell equations in vacuum can be decomposed into plane waves 
with positive and negative frequencies. This decomposition gives the following Fourier 
representation of F(r,t): 

F(r, t) = VN j e(fc) [f L (k) e -^ t+ik - r + f^e^-^} , (16) 

where the complex polarization vector e(k) — [li(k) + il2(k)}/^/2 has the following 
properties: 

ck x e(k) = — icjfc e(fc), (17a) 
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e(fc)-e(fc) = 0, (176) 

e*(k)-e(k) = 1, (17c) 

e*(k) x e(k) = in k , (17 d) 

e*(k)-e(-k) = 0, (17e) 

e(k)xe(k) = 0, (17 f) 

e*(k) ej (k) = i (Sij + J ■ {VI g) 



The identification of the Fourier coefficients with the components of the photon wave 
function in the formula (|16|) will be justified in the next section where we will unify 
the field picture and the photon picture. The second term in (I16p involves complex 
conjugation. This is dictated by the fact that the photon energy is always positive. 
Therefore, the time evolution of the wave function is given by the factor exp(— iw^t). 
Therefore, the reversal of the sign in the exponent requires complex conjugation. We 
pulled out the factor y/N to assure the normalization of f. 



4. Separation of angular momentum 

We shall combine now the field picture and the photon picture to obtain the 
decomposition of the total angular momentum of the field. To this end, we substitute 
the Fourier representation of the field into the formulas (|15al) - (jl5ci|) . 

H = N ^f(k)-huj k f(k), (18a) 



k 



d 3 k 

P = N —f(k)-hkf(k), (18b) 



J = N f P^f(k)-[ihDxk + hxn k ]f(k), (18c) 
J nio k 

/d 3 k 

—f(k)-ihu: k Df(k). (18d) 
nw k 

Note, that the resulting expressions have the form of quantum mechanical expectation 
values 

H = N(f\H\f), (19a) 

P = N(f\P\f), (19b) 

J = N(f\J\f), (19c) 

K = N(f\K\f). (19d) 

These formulas exhibit a perfect agreement between the results obtained from the 
particle picture and from the field picture, as Darwin had anticipated. Every value 
calculated for the total electromagnetic field is a product of of the quantum mechanical 
average value per one photon, multiplied by N. That means that the normalization 
factor N is the total number of photons. We may now unambiguously split the total 
angular momentum of the electromagnetic field (|18cp into two parts. The vector J a 
whose integrand is perpendicular to the wave vector is the orbital part and the vector 
J s with integrand parallel to the wave vector is the spin part represented by helicity, 

/d 3 k 
—fUk)-[ihDxk]f(k), (20a) 
huj k 
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Js=N 



/d 3 k f d 3 k 

—f(k)-hxn k f(k) =N J —n k [\f L {k)\ 2 - \f R {k)\ 2 } , (206) 

The final step of our analysis is the proof that the expressions for J Q and J s coincide 
with those obtained by Darwin. To this end, we employ the relation between E(k) 
and f(fe) that follows from the formulas (|3J and ([To]) 

E(k) = ^ [e(k)f L (k) + e*{k)f R {k)] . (21) 

Upon substituting this relation into the second term in ([2]), with the use of the 
properties of the polarization vectors (|17ri|) and (17/1, we obtain 

/d 3 k 
— E*(k) x E(k) 

= - iN J^n t e *( fc ) x < k )\fL{k)\ 2 + e(k) x e*(k)\f R (k)\ 2 ] = J s . (22) 

In the same way we may establish the equality of the orbital part in the Darwin form 
and in quantum mechanics of photons. Note that the separation of the total angular 
momentum into its orbital and spin parts is conserved in time since both parts are 
separately time independent. 

As an illustration, we consider the Bessel beam characterized by the frequency 
c\k\, the z-component of the total angular momentum m, the component k z of the 
wave vector in the z-direction, and the helicity ±1. In this case, the Darwin vector 
(|21[) (up to a normalization factor) as given in [33J has the form 

(— (fc 2 /fc) cos0 ± isin0 \ 
-(fc 2 /fc)sin(/)Ticos0 e im ' t 'S(k^-k , ± )5(k z -k' z ).(23) 
{kjkf ) 

Since the Bessel beam has an infinite extension in space, both parts of the total angular 
momentum are infinite. However, their ratio is well defined. For the components in 
the beam direction, the ratio of the orbital to spin parts equals to mk/k z =p 1. 

In order to express J and J s as integrals in coordinate space we have to invert 
the Fourier transformation in (J3j for t = as follows: 

d 3 r 



E(r) + x B(r) 



(24) 



2(2tt) 3 / 2 ' 

where we made use of Maxwell equations. Inserting this formula into the Darwin 
expression for the spin part J s , after the integration over fc, we obtain 

where we used the formula 
d 3 k e lk " 



(2tt) 3 |fe| 2 iw\r\' ^ 
This gauge-invariant integral representation of J s becomes equal to the last term in 
(JXJ) if the vector potential is identified with the following integral: 

_ VxBM. (27) 

J 47r \r - r'\ 

This representation of the vector potential is valid in the transverse gauge, as has 
been anticipated. Note that the seemingly local form of the formula ([T]) is misleading 
because the gauge invariant vector potential is a nonlocal function of the magnetic 
field. 
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5. Conclusions 

We have shown that the separation of the total angular momentum of the 
electromagnetic held into its orbital and spin parts dictated by quantum mechanics 
of photons reproduces the results derived from the properties of Maxwell fields by 
Darwin. This separation, when expressed in the form of coordinate-space integrals, 
coincides with the results derived heuristically by many authors, provided the vector 
potential is related to the magnetic held by the integral formula ([27)1 . In contrast to 
energy, momentum, and the total angular momentum of the electromagnetic field, the 
orbital angular momentum and the spin parts cannot be expressed as integrals of local 
densities: they are intrinsically nonlocal objects. 
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